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Using a viscosity-based survival scale for geometrical perturbations formed in the early stages of 
relativistic heavy-ion collisions, we model the radial flow velocity during freeze-out. Subsequently, we 
employ the Cooper-Frye freeze-out prescription, with first-order viscous corrections to the distribu¬ 
tion function, to obtain the transverse momentum distribution of particle yields and flow harmonics. 

For initial eccentricities, we use the results of Monte Carlo Glauber model. We fix the blast-wave 
model parameters by fitting the transverse momentum spectra of identified particles at the Large 
Hadron Collider (LHC) and demonstrate that this leads to a fairly good agreement with transverse 
momentum distribution of elliptic and triangular flow for various centralities. Within this viscous 
blast-wave model, we estimate the shear viscosity to entropy density ratio r]/s ~ 0.24 at the LHC. 
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I. INTRODUCTION 

High-energy heavy-ion collision experiments at the 
Relativistic Heavy Ion Collider (RHIC) [1, 2] and the 
Large Hadron Collider (LHC) [3-5] have conclusively es¬ 
tablished the formation of a strongly interacting Quark- 
Gluon Plasma (QGP). The QGP formed in these colli¬ 
sions exhibit strong collective behaviour, and therefore 
can be studied within the framework of relativistic hy¬ 
drodynamics. The hydrodynamical analyses of the flow 
data suggests that the QGP behaves like a nearly perfect 
fluid with an extremely small shear viscosity to entropy 
density ratio rj/s [ 6 , 7]. Local pressure gradients due to 
deformations and inhomogeneities in the the initial stages 
of the collision results in anisotropic fluid velocity. These 
anisotropies subsequently translate into flow harmonics 
describing the momentum asymmetry of produced par¬ 
ticles [ 8 , 9]. The 77 /s of the fluid governs the conversion 
efficiency, which results in a suppression of elliptic flow 
and higher order flow coefficients [10-23]. 

Apart from hydrodynamics, the observables pertaining 
to collective behaviour of QGP can also be studied using 
the so-called blast-wave model. Using a simple functional 
form for the phase-space density at kinetic freezeout, 
Schnedermann et. al. [24] approximated hydrodynamical 
results with boost-invariant longitudinal flow [25]. They 
used this blast-wave model to successfully fit the trans¬ 
verse momentum spectra with only two parameters: a ki¬ 
netic temperature, and a radial flow strength. However, 
this model was only valid for central collisions at midra¬ 
pidity. In order to make it applicable for non-central 
collisions, Huovinen et. al. [26] generalized this model to 
account for the anisotropies in the transverse flow pro¬ 
file by introducing an additional parameter. This new 
parameter controlled the difference between the strength 
of the flow in and out of the reaction plane. This lead 
to a fairly good 6 t with the measured elliptical flow as a 
function of transverse momentum. However, the STAR 


Gollaboration achieved better fits when they generalized 
the model even further by introducing a fourth param¬ 
eter to account for the anisotropic shape of the source 
in coordinate space [27]. Teaney made the first attempt 
to estimate the effect of viscosity on elliptic flow using 
a variant of the blast-wave model model [28]. However, 
the centrality dependence of the ht parameters left the 
model with little predictive power. 

In this paper we generalize the blast-wave model to in¬ 
clude viscous effects by employing a viscosity-based sur¬ 
vival scale for geometrical anisotropies, formed in the 
early stages of relativistic heavy-ion collisions, in the 
parametrization of the radial flow velocity. The present 
model has hve parameters, including 77 /s, which has to be 
fitted only for one centrality. In the Gooper-Frye freeze¬ 
out prescription for particle production, we consider the 
first-order viscous corrections to the distribution function 
[28]. In essence, we provide a model which incorporates 
the important features of viscous hydrodynamic evolu¬ 
tion but does not require to do the actual evolution. We 
use this viscous blast-wave model to obtain the trans¬ 
verse momentum distribution of particle yields and flow 
harmonics for LHC. The blast-wave model parameters 
are fixed by fitting the transverse momentum spectra of 
identihed particles. Subsequently, we show that this leads 
to fairly good agreement with transverse momentum dis¬ 
tribution of elliptic and triangular flow for various cen¬ 
tralities as well as centrality distribution for integrated 
flow. We estimate the shear viscosity to entropy density 
ratio 77 /s ~ 0.24 at the LHC, within the present model. 


II. BLAST WAVE MODEL 

The blast-wave model has been used extensively to fit 
experimental data and it provides good description of 
spectra and elliptic flow observed in relativistic heavy-ion 
collisions [26-30]. The previously used blast-wave mod- 
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els employ a simple parametrization for the flow velocity 
of boost invariant ideal hydrodynamics. The most im¬ 
portant feature is the parametrization of the transverse 
velocity which is assumed to increase linearly with re¬ 
spect to the radius. This parametrization is found to be 
in agreement with hydro results [31, 32]. This essentially 
leads to an exponential expansion of the fireball in the 
transverse direction, hence the term blast-wave. Apart 
from boost invariance, the model also assumes rotational 
invariance. In the following, we quickly outline the key 
features of the blast-wave model. 

In order to consider a boost invariant framework, it is 
easier to work in the Milne co-ordinate system where. 


T ^ \/P - 

(1) 

rjs = tanh“^(z/t), 

(2) 

r = \Jx'^ -\- y'^, 

(3) 

If = atan 2 ( 2 /, x). 

(4) 


The metric tensor for this co-ordinate system is = 
diag{l, —T^, —1, —r^). Boost invariance and rotational 
invariance implies = 0, whereas linearly rising 

transverse velocity flow profile leads to vT ^ r. The blast- 
wave model further assumes that the particle freeze-out 
happens at a proper time r/ having a constant temper¬ 
ature Tf and uniform matter distribution, in the trans¬ 
verse plane. In summary, the hydrodynamic fields are 
parametrized as [28] 


T = TfQ[R-r), 

(5) 

= uo^ 0(i? - r) 

H 

(6) 

= 0, 

(7) 

u" = Jl + 

(8) 


where R is the transverse radius of the fireball at freeze- 
out. The expression for u'^ is obtained by requiring that 
the fluid four-velocity satisfy the condition = 1. 

The hadron spectra can be obtained using the Cooper- 
Frye prescription for particle production [33] 

(9) 

where is the oriented freeze-out hyper-surface and 
f{x,p) is the phase-space distribution function of the 
particles at freeze-out. The distribution function can be 
written in terms of the equilibrium and non-equilibrium 
parts, / = fo + Sf. The equilibrium distribution function 
is given by 

exp{u^pf^/T) + a ’ 

where a = -1-1 for baryons and a = — 1 mesons. 

For small deviations from equilibrium, i.e., Sf ^ /o, 
we use the Grad’s 14-moment approximation for the non¬ 
equilibrium part [34, 35] 


where fo = l — afo and tTq,^ is the shear stress tensor. 
Approximating the shear stress tensor with its first-order 
relativistic Navier-Stokes expression, tTq,^ = 2rf\/(^^.u, 
the expression for the 14-moment approximation reduces 
to 


= ( 12 ) 

Here rj is the coefficient of shear viscosity, s = {e + P)/T 
is the entropy density and the angular brackets denote 
traceless symmetric projection orthogonal to the fluid 
four-velocity [36]. The form of in the case 

of blast-wave model is calculated in Appendix A. 

The anisotropic flow is defined as 


dcf) cos[n{(j) — d^n)] 


dN 


VniPr) = 


dypt dpT d4> 


d(f> 


dN 


(13) 


dypt dpT d4> 


where d'n is the n-th harmonic event-plane angle. In the 
present case, we do not consider event-by-event fluctua¬ 
tions and therefore = 0. Up to first order in viscosity 

|28|, 


1 - h 1 


d(j) ■ 


dypt dpT dfj) 

dN^^^ 

d(j) cos[n{(j) — 4'„)] 




-k 


dypt dpT d(j) 


d4> 




(14) 


dypt dpT d4> 


where the superscript ‘(0)’ denotes quantities calculated 
using the ideal distribution function, Eq. (10), and ‘(1)’ 
denotes those obtained using the first-order viscous cor¬ 
rection, Eq. (12). 


III. VISCOUS BLAST-WAVE MODEL 


The definition of the participant anisotropies, via 
the Fourier expansion for a single-particle distribution is 


/(¥>) 


1 

27T 


oo 

1 + 2^ £„ cos[n((^- V'n)] , 

n—1 


(15) 


where V'n are the angles between the x axis and the ma¬ 
jor axis of the participant distribution. The geometrical 
anisotropies in the initial particle distribution, £„, even¬ 
tually converts to anisotropies in the radial fluid velocity. 


r 



OO 

1 + 2 ^ cos[n(v? - -ipn)] 

n—1 


( 16 ) 
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FIG. 1: The harmonics form standing waves on the fireball 
circumference having radius R at freezeout. 


In the following, we determine the conversion efficiency 
of the initial eccentricity to anisotropy in the radial fluid 
velocity, u„/£„. 

Using the well known dispersion relation for sound in 
a viscous medium [35] , 

uj = Csk + ik'^^ , (17) 


the authors of Ref. [37] introduced a viscosity-based sur¬ 
vival scale which all structures formed by point like per¬ 
turbations should attain at freeze-out. In the above equa¬ 
tion, r] is the coefficients of shear viscosity and Cs is the 
speed of sound in the medium. In the present work, we 
ignore the contribution due to bulk viscosity. Using a 
plane-wave Fourier ansatz, exp(*a;< — ikx), we observe 
that the amplitudes of the stress tensor harmonics with 
momentum k are attenuated by a factor 


k) = exp 



kH' 
~Y. 


STf^’^{0,k), 


(18) 


where we have suppressed the oscillatory pre-factor. We 
note that the presence of momentum squared in the ex¬ 
ponent leads to enhanced effect of viscosity for the higher 
harmonics. We expect the same qualitative behaviour for 
the radial flow velocity as will be explained in the follow¬ 
ing. 

First and foremost, we notice that each harmonics 
is essentially a damped oscillator with wave-vector k. 
Moreover, throughout the evolution the harmonics form 
standing waves on the fireball circumference, as shown 
in Fig. 1, whose amplitude is progressively damped due 
the viscous effects. Therefore, the fireball circumference 
is an integer multiple of the wavelength with wave-vector 
k, i.e.. 


2ttR = n —, (19) 

k 

where R is the transverse radius of the fireball at freeze¬ 
out. Hence, at the freeze-out time tf, the wave amplitude 
reaction is given by 




= exp 




R^Tf 


( 20 ) 


where T/ is the freeze-out temperature. 

In absence of viscosity, the initial geometrical pertur¬ 
bations in the fluid will result in the development of ra¬ 
dial flow velocity and the conversion efficiency will re¬ 
main the same for all harmonics. In the case of a viscous 
medium, however, the conversion efficiency of the initial 
geometrical perturbation to radial fluid velocity must be 
proportional to the wave amplitude reaction. 


— = ao exp 




R^Tf 


( 21 ) 


where ao is the constant of proportionality. The sudden 
stopping of the damped oscillator at the freeze-out time 
may lead to certain phases due to the oscillatory pre¬ 
factor. These phases can, in general, lead to secondary 
peaks in the power spectrum of higher harmonics. How¬ 
ever, as no secondary peaks has been observed in the 
spectrum of relativistic heavy-ion collisions, we will con¬ 
tinue to ignore these phase factors. 

We emphasize that the acoustic damping should be 
applied to the hydrodynamic variables, such as the mo¬ 
ments of the flow velocity, rather then to the final 
state observables such as as done in Ref. [38]. In 
Ref. [38], Shuryak and Zahed (S-Z) proposed that the ra¬ 
tio of the initial eccentricity e„ to the final pr-integrated 
Vn should be proportional to the wave amplitude reac¬ 
tion, i.e., the r.h.s. of Eq. (21) should be equal to u„/£„. 
However, £„ is the eccentricity in the configuration space 
whereas is the momentum anisotropy of the parti¬ 
cles after freeze-out. The momentum space asymmetries 
are not hydrodynamic variables and are only affected in¬ 
directly via damping. On the other hand, the acous¬ 
tic damping should be applicable to hydrodynamic vari¬ 
ables and it should only capture the viscous effects of the 
hydrodynamic evolution. This assumption also misses 
the additional effect of viscosity at freeze-out using the 
Cooper-Frye formula. Moreover, it does not provide the 
opportunity to study the pt dependence of anisotropic 
flow and hence to estimate the viscosity of the expand¬ 
ing medium. 


IV. INITIAL CONDITIONS 

In this section, we set-up the initial conditions of the 
collisions in order to reduce the number of free param¬ 
eters in the blast-wave model. To this end, we evaluate 
the parameters corresponding to the initial geometry of 
the collisions. We also approximate the subsequent trans¬ 
verse expansion of the fireball by using the radial velocity 
parametrization in the blast-wave model. 

We consider the collision of two identical nuclei with 
mass number A. The radius of each nucleus is given 
by i?o = 1.25A^/^ fm and the impact parameter is b, 
as shown in Fig. 2. The shaded region in Fig. 2 is the 
overlap zone of the colliding nuclei. We draw a circle 
of radius Tq, with its centre coinciding with that of the 
overlap zone, in such a way that the boundary of the 
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FIG. 2: Initial transverse geometry of the collision of two 
identical nuclei with radius Rq and impact parameter b. The 
shaded region represents the the overlap zone of the colliding 
nuclei. The circle with radius ro is drawn such that it equally 
divides the boundary of the overlap zone in four parts, i.e., 
123 = 345 = 567 = 781. 



FIG. 3: (Golor online) Gentrality dependence of the initial 
energy density Ci and freeze-out time r/ for Pb-|-Pb collisions 
at y/SNN = 2.76 TeV scaled by the corresponding values in 
0 — 5% central collisions. The inset shows centrality depen¬ 
dence of the freeze-out radius of the fireball R, obtained using 
Eq. (24). 


overlap zone is equally divided in four parts, i.e., 123 = 
345 = 567 = 781; see Fig. 2. This approximates £2 as 
the second harmonics of initial geometrical fluctuations, 
analogous to the n = 2 case as shown in Fig. 1. The 
radius tq is therefore the initial transverse radius of the 
expanding fireball and is given by 

ro = Ub^-2bRo^j2 + ^^+4R^^ , ( 22 ) 


For the ideal hydrodynamic evolution of relativistic 
fluid, in the one-dimensional boost-invariant scenario, the 
evolution of the energy density follows e oc Assum¬ 

ing the initial thermalization time and the final freeze-out 
energy density (i.e., the freeze-out temperature) to be 
same for all collisions, we get 

rf=rfo[^) ■ (25) 


which reduces to rg = Rq for head-on collisions (6 = 
0). Since £2 is the most prominent eccentricity for non¬ 
central collisions, all other geometrical eccentricities are 
treated as boundary perturbations to this circle. 

The subsequent transverse expansion of the fireball is 
obtained by employing the radial velocity parametriza- 
tion of the blast-wave model. Using the perturbation-free 
expression for the transverse velocity, Eq. (6), we get 


^ dr r 



(23) 


After performing the straightforward integration, we ob¬ 
tain a transcendental equation for the freeze-out radius, 

R, 


i? = roexp(^^^^) , (24) 

which can be solved for R given the isotropic expansion 
velocity ug and the freeze-out time tq. In the follow¬ 
ing, using Bjorken’s scaling solution for one-dimensional 
boost-invariant expansion, we obtain an expression to 
determine the freeze-out times for non-central collisions 
once it is fixed for the central one. 


Here r/o is the freeze-out time for most central collisions, 
which has to be fixed by fitting the corresponding trans¬ 
verse momentum spectra. The freeze-out times for other 
centralities can then be obtained using the above equa¬ 
tion and therefore they are not free parameters. The 
ratio (-i/tig is the initial central energy density scaled by 
its corresponding value in most central collisions. 

Figure 3 shows tiftig for various centralities obtained 
using the Glauber model calculations [39, 40] in the case 
of Pb-Pb collisions at = 2.76 TeV at the LHC. 

We observe that tijtig (and hence the initial temper¬ 
ature) decreases for non-central collisions compared to 
central ones. Therefore according to Eq. (25), freeze-out 
happens earlier in peripheral collisions which is also re¬ 
flected in Fig. 3 for Tfjrfg. The inset of Fig. 3 shows the 
centrality dependence of the freeze-out radius of the fire¬ 
ball R, obtained using Eq. (24). We find a rather large 
transverse radius of the fireball at freeze-out. Finally, the 
parameters that we need to fix within the viscous-blast 
wave model to fit the spectra are the freeze-out temper¬ 
ature Tf, the freeze-out time for central collision Tfg and 
the unperturbed maximum radial flow velocity ug. An in¬ 
terplay of the coefficient of proportionality ag in Eq. (21) 
and r]/s will be crucial to reproduce the flow harmonics. 
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FIG. 4: (Color online) Transverse momentum distribution of 
particle multiplicities in Pb+Pb collisions at = 2.76 

TeV. We show results for tt"*", K'^ , and p in two centrality 
ranges, (a): 0 — 5% and (b): 20 — 30%. The symbols represent 
ALICE data [41] at midrapidity and the lines correspond to 
viscous blast-wave calculations. 


V. RESULTS AND DISCUSSIONS 

In this section, we show our results for Pb-|-Pb colli¬ 
sions at y/SNN = 2.76 TeV and compare it with experi¬ 
mental data measured at LHC by the ALICE and ATLAS 
collaborations. 

Figure 4 shows the transverse momentum distribution 
of pions, kaons, and protons spectra for 0 — 5% and 
20 — 30% central Pb-t-Pb collisions at ^/SNN = 2.76 TeV 
measured at LHC by the ALICE collaboration (symbols) 
at midrapidity [41] and calculated using the viscous blast- 
wave model (lines). The results are obtained using root 
mean square values of eccentricities £2 and £3 in a Monte- 
Carlo Glauber model, with a shear viscosity to entropy 
density ratio = 0.24. We observe that the spectra 
for TT'*' and K'^ from the viscous blast-wave model are in 
good overall agreement with the experimental data for a 
freeze-out temperature of 120 MeV. On the other hand, 
within the viscous blast-wave model, the proton yield for 
a freeze-out temperature of 120 MeV is severely under¬ 
estimated. To obtain an overall fair agreement with the 
experimental data, the freeze-out temperature for pro¬ 
tons is considered to be 135 MeV. The freeze-out time 
for 0 — 5% most central collision was found to be 8 fm. 
For other centralities, the freeze-out time was obtained 
by using Eq. (25) and results shown in Fig. 3. 

Figure 5 shows our results for the Vnipr), in compar¬ 
ison with the ATLAS data [4] for various centralities. 
We find overall fair agreement with the data for elliptic 
(^ 2 ) and triangular (v^) flow, at all centralities. This is 
achieved by choosing a single fixed value rj/s = 0.24 to 
obtain the required suppression (relative to ideal blast- 
wave results) of V 2 and U 3 , for all centralities. Apart from 
ry/s, flow also depends on the constant of proportional¬ 
ity ao appearing in Eq. (21). It controls the conversion 



FIG. 5: (Color online) Transverse momentum dependence of 
the anisotropic flow coefficients Vn(pT) of charged hadrons, for 
n = 2 and 3, calculated at various centralities in Pb-|-Pb col¬ 
lisions at ^ysNN = 2.76 TeV in the viscous blast-wave model 
(lines) with pjs = 0.24 as compared to the ATLAS data [4] 
(symbols). 


efficiency of the initial eccentricity to final fluid veloc¬ 
ity. We consider the initial eccentricity £„ to be the rms 
values of the eccentricities obtained in the MC-Glauber 
model, as given in Ref. [42]. A large value of Oq means 
larger conversion of eccentricity leading to increased flow 
velocity and hence higher V 2 and V 3 . On the other hand, 
as is well known, an increase in ry/s leads to suppression 
of V 2 and V 3 . Large (small) value of Oq can be compen¬ 
sated by choosing a higher (lower) value of ?y/s up to a 
certain extent. However, beyond a certain range of ry/s, 
the relative behaviour of V 2 and V 3 is destroyed. 

In order to match the elliptic and triangular flow data, 
we find the most suitable parameter values for ao = 0.4 
and rjjs = 0.24. Moreover, we find that Vn is insensi¬ 
tive to a for transverse velocity of the form ~ (rIK)°‘. 
This may be attributed to the fact that the exponent a 
controls the rate of isotropic transverse expansion. The 
anisotropic flow originates from the initial eccentricity 
which translates into final flow. Therefore is sensi¬ 
tive to Un which depends on ag and ?y/s, as is apparent 
from Eq. (21). On the other hand, it should be noted 
that the slope of the particle spectra is sensitive to the 
exponent a and could be tuned to get a better ht with 
the experimental data. However, in the present work, we 
are interested in the anisotropic flow and therefore we set 
0 = 1 and do not fit it to match the particle spectra. 

In Fig. 6(a), pT-integrated values of V 2 and V 3 obtained 
from the viscous blast-wave model are compared with 
ALICE data [3], as a function of centrality. We observe 
that using the same constant rj/s = 0.24 for all central- 
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Centrality (%) 


FIG. 6: (Color online) (a): Centrality dependence of the pr 
integrated anisotropic flow coefficients of charged hadrons 
in Pb+Pb collisions at ^/snn = 2.76 TeV calculated in the 
viscous blast-wave model (lines) with rj/s = 0.24, as compared 
to ALICE data [3] (symbols), (b): Centrality dependence of 
the ratio Vnlsn in the viscous blast-wave model where £n 
is the rms values of the eccentricities obtained in the MC- 
Glauber model [42]. In panels (a) and (b), we show results 
for n = 2 and 3. (c): Centrality dependence of the ratio 

U 2 /U 3 for ALICE data (symbols), using the viscous blast- 
wave model (solid line) and from the Shuryak-Zahed estimate 
(dashed line). 


ities, the model shows a good agreement with the data. 
Since V 2 is driven mostly by the initial spatial anisotropy, 
it exhibits a strong centrality dependence compared to ^ 3 . 
Figure 6 (b) shows the conversion efficiency of the initial 
spatial anisotropy into the final momentum anisotropy. 
A linear relation between Vn and as observed in some 
hydrodynamic calculations [23], is not obtained within 
the viscous blast-wave model presented here. On the 
other hand, in view of the non-linear nature of the hy¬ 
drodynamic equations a linear relation between Vn and 
Sn is not obvious. Indeed other calculations [43] as well 
as a recent analysis of the LHC data [44] also result in 
a similar centrality dependence as obtained in our anal¬ 
ysis In Fig. 6 (c), we show the ratio V 2 /V 3 as a function 
of centrality for the ALICE data (symbols), the present 
viscous blast-wave model (blue solid line) and the esti¬ 
mate due to Shuryak and Zahed (red dashed line). We 
see that the viscous blast-wave model provides a better 
agreement with ALICE data compared to the S-Z esti¬ 
mate. However, we should keep in mind that the freeze- 
out parameters used in the S-Z estimate is the same as 
that of the viscous blast-wave fit values. 


VI. CONCLUSIONS AND OUTLOOK 

In this paper we have generalized the blast-wave model 
to include viscous effects by employing a viscosity-based 
survival scale for geometrical anisotropies formed in the 
early stages of relativistic heavy-ion collisions. This vis¬ 
cous damping is introduced in the parametrization of 
the radial flow velocity. The viscous blast wave model 
presented here involved five parameters, including 77 /s, 
which has to be fitted for only one centrality. This 
model therefore incorporates the important features of 
viscous hydrodynamic evolution but does not require to 
do the actual evolution. We have used this viscous blast- 
wave model to obtain the transverse momentum distri¬ 
bution of particle yields and anisotropic flow harmonics 
for LHC. The blast-wave model parameters were fixed 
by fitting the transverse momentum spectra of identified 
particles. We demonstrated that a fairly good agreement 
was achieved for transverse momentum distribution of el¬ 
liptic and triangular flow for various centralities as well 
as centrality distribution for the integrated flow. Within 
the present model, we estimated the shear viscosity to 
entropy density ratio rj/s ~ 0.24 at the LHC. 

One of the drawbacks of the present model is that 
we have employed root mean squared eccentricity over 
number of events, which is analogous to “single shot” 
hydrodynamic evolution. On the plus side, the present 
model could also be implemented on an event-by-event 
basis. Another difficulty that we encountered was that 
in order to fit the proton spectra, we had to consider a 
different freeze-out temperature (135 MeV) compared to 
that for pions and Kaons (120 MeV). This problem could 
be addressed by parametrizing the transverse velocity in 
the form u” ~ (r/i?)“ and fitting a for different particle 
species, separately. We leave these problems for a future 
work. 
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Appendix A: Viscous stress tensor 

In this Appendix, we calculate the viscous tensor 
and therefore obtain the viscous corrections to 
the distribution function at freeze-out. We work in Milne 
co-ordinate system, Eqs. (l)-(4) with the metric ten¬ 
sor ( 7 ^ 1 , = diag(l, —r^, —1, —r^). Therefore, the inverse 
metric tensor is = diag(l, — 1 /t^, — 1 , —1/r^), its de¬ 
terminant g is y/—g = Tr and the non-vanishing Christof- 
fel symbols are F”^^^ = r, F^j^^ = 1/r, F”,^ = —r, and 
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= l/r. Using the parametrization of the fluid veloc- where we have used the fact that = 1 + 

ity given in Eqs. (6)-(8), we get 


A’'‘^ = 0, A‘^^ = —A" = -l-(u'■)^ (Al) 


d^vJ' = —, =-^^^^nu„sin[n((/?-■!/'«)]• 


n—1 


(A2) 


where A^"^ = — u^u'' is the projection operator or¬ 

thogonal to the fluid four-velocity. 

To fix the time derivatives of the fluid velocity, we as¬ 
sume that if the particles are freezing-out, they are free 
streaming, which means that Du^ = 0. Here D = 
is the co-moving derivative and is the covariant deriva¬ 
tive. With this prescription, we have 


dru'f’ = 0, 

= —V drVj' = — 


r^2 


drVd' = V drVd’ = — 


ru^ 

{v^ 

r{u^Y 


(A3) 


where v = vT/v7 is the radial velocity. The expansion 
scalar is given by 

1 _ 7 /^ nj^ 

=-1-h + drU^ + drU", 

V-9 T r ^ 

u’^ (un^ ,,,, 

—-1" 2- , .2' 

r r ryu^y 

Assuming boost invariance, the spatial components of the 
viscous tensor are given by 

- -u^Du^ - -u‘^Du^ 


CXJ 

= Sin[n((p - y ^)], 


(AS) 


n—1 


= _ a yV- yu'^Du'f’ 

r 


3 V-9 


1 




3 It r r{u^y] 


(A6) 


yir^r) ^ _ 1 ^rr 

3 V=5 




— + 
r r 


(A7) 


3 y/—g 

2 ^-2-- 


3 L r r r(w'^) 


^{ry^Vs) ^ ^{‘PyVs) ^ 0 . 


(AS) 


(A9) 


To obtain the temporal components of the viscous 
stress energy tensor, we use the Landau frame condition, 
= 0 . 

V^^U^Vr+V^^u''Vr=0 - 


(AlO) 

yiryVs) = 0, (All) 
(A12) 
(A13) 

Therefore, from Eqs. (AlO) and (A12), we see that 




{uV'^u^ {uV^ (uO® 


3 L T 


r{u‘^y r{u‘^y. 


(A14) 


Next, in order to verify our algebra, we confirm that the 
viscous stress tensor is traceless, i.e., = 0. 

Using Eqs. (A6), (A7), (A8) and (A14) 


[vJ'f [vdf 


3 L 
1 

~ 3 
1 


r{u^y r{u^y. 


2 ^- 2 --^ 
r T r{u^)‘^. 

r 


3 L r 


1 


u u 


(un 


r 

r\3 


3 It r r{u^Y\ 

= 0. (A15) 

For a particle at the space-time point (t, 77 s, r, tp) 
with the four momentum = (E, , p^) = 

(rriT coshy, pt cos Pp, px sin(/?p, mr sinhy), we get 

p^ = rriT cosh{y - gs) ^ Pr = wt cosh(y - 77 ^), 

rpVs _ sinh(7/ - 77 ^) ^ Pr,, = -t mr sinh(7/ - Ps), 

p^ = Pt cos(pp -p) ^ pr = -pT cos(pp - p), 

r p'^ = Pt sin(pp - p) => p^ = -rpr sin(pp - p). 

(A16) 

The oriented freeze-out hyper-surface is dSp = 
(Tdps fdr dp, 0 , 0 , 0 ), and therefore the integration mea¬ 
sure is given by 

p^dEp = rriT cosh(y — 77 ^) Tdps rdr dp. (A17) 









































The viscous correction to the equilibrium distribution 
function is proportional to 

+ pIV^^u^'^ + 2p,prV^^u^'> + 2prP^V^^u^'> 

+ 2prP^V^^u‘^K (A18) 

The final form of / = fo+Sf obtained from Eq. (12) using 


the above equation is required to evaluate the spectra 
given by 
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(A19) 
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